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^ . Abstract 

o ■ 

I We consider Brans-Dicke theory with a self-interacting potential in Einstein conformal 

frame. We show that an accelerating expansion is possible in a spatially flat universe for 
large values of the Brans-Dicke parameter consistent with local gravity experiments. 
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Recent observations on expansion history of the universe indicate that the universe is experi- 
encing a phase of accelerated expansion [1]. There are two classes of models in the literature 
I for explaining these observations. In one class, cosmic acceleration is attributed to some exotic 
matter components with a negative pressure, dubbed dark energy [2] . This class also includes 
quintessence models [3], a minimally coupled scalar field with a potential which can generate 
a sufficiently negative pressure at the present epoch. It should be pointed out that in an ob- 
servational point of view, models with a cosmological term or adiabatic ACDM (cosmological 
k>( ; constant or non-evolving dark energy) seem to be in agreement with cosmological observations 
' [4]. However, the observational upper bound differs from theoretical expectations by more than 
100 orders of magnitude [5] . This large discrepancy avoids people to consider the cosmological 
constant viable candidate of dark energy [6] . 

Alternatively, there is a class of models that propose a modification of the gravitational part 
of the Einstein-Hilbert action. This includes scalar-tensor theories, scalar fields non-minimally 
coupled to gravity. Here we shall consider a self-interacting Brans-Dicke (BD) theory [7] as a 
prototype of scalar-tensor theories. The original motivation of the BD theory was the search 
for a theory containing Mach's principle which has found a limited expression in general rel- 
ativity. As the simplest and best-studied generalization of general relativity, it is natural 
to think about the BD scalar field as a possible candidate for producing cosmic acceleration 
without invoking a quintessence field or exotic matter systems. In fact, there have been many 
attempts to show that BD model can potentially explain the cosmic acceleration. It is shown 
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that this theory can actually produce a non-decelerating expansion for low negative values of 
the BD parameter uj [8]. This conflicts with the lower bound imposed on this parameter by 
solar system experiments [9]. Some authors propose modifications of the BD model such as 
introducing a potential function for the scalar field [10], or considering a field-dependent BD 
parameter [11], without resolving this problem. 

All the works in this context use Jordan frame representation of BD theory. It is however well- 
known that this theory, like any other scalar-tensor theories, can be represented in the so-called 
Einstein frame by using a conformal transformation [12] [13]. Although these two conformal 
frames are mathematically equivalent there are some debates on their physical equivalence. 
Here there is a point which should be made clear. Despite the fact that the problem of phys- 
ical status of the two conformal frames is open there is a tendency in the literature to ignore 
this problem and to work in Jordan conformal frame. The reason may be related to reluctance 
in accepting of the violation of weak equivalence principle due to anomalous coupling of the 
scalar field to matter systems in Einstein frame. It is however important to keep in mind 
that the physical metric should be singled out already in the vacuum sector of the theory and 
the coupling of a given metric to matter systems is determined by the physical significance 
ascribed to it [14]. Thus a criterion based on the coupling of matter with gravity would be 
effective only if the physical frame were determined on an independent ground. Apart from 
this point, anomalous gravitational coupling in Einstein conformal frame does not necessarily 
mean violation of weak equivalence principle. There is still a possibility that the effective mass 
of the scalar field be scale dependent. In this chameleon mechanism [15], the scalar field may 
acquire a large effective mass in Solar System scale so that it hides local experiments while at 
cosmological scales it is effectively light and can provide an appropriate cosmological behavior. 
Along these lines, we would like to consider Einstein frame formulation of the theory as a 
representation which provides different possibilities in a cosmological setting with respect to 
the Jordan frame. In particular, we will focus on the question that whether it is possible to 
achieve accelerating expansion of the universe for sufficiently large values of the parameter u. 
We will show that the answer is affirmative within a class of solutions which corresponds to a 
specific form of the potential function of the BD scalar field. 

We begin with a modified form of the BD action in Jordan frame 

SjF = j d^x^m - ^^'^'^V^^V.^ - y (0)) + Sm{g^,., ^) (1) 

where (p is the self-interacting BD scalar field with a potential function V{(f)), w is a constant 
parameter and is the action of matter which depends on the metric g^^, and some matter 
fields collectively denoted by -0. A conformal transformation 

Qij-u Qfiu = ^'^Qu.i. (2) 

with ft — brings the above action into the Einstein frame [12] [13]. Then a scalar field 

redefinition 
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with 00 ~ C ^, > and a; > — | transforms the kinetic term of the scalar field into a 
canonical form. In terms of the variables {g/xu, ip) the BD action in the Einstein frame is [12] 
[13] 

Sef = / d'xV^ij^ - Ig'^-V.^V.^ - U{ip)) + Srn{g,,, ^) (4) 

where 

Sm = j d'^xy/^ exp(-8J ^^^ y?) Lrn{gtj,u,ij) (5) 
Here V ^ is the covariant derivative of the rescaled metric g^j^v and 



U{^) = V{<P{ip)) exp(-8y ^^V) (6) 
is the Einstein frame potential. 

Variation of the action (4) with respect to gfj,i, and (p leads to the following field equations 

G,, = 87rGiT,, + Tp (7) 

^ d(p 2 ^ ^ 

where 

= V^(^V,(^ - ^^^.V^(^V> - U{ip)g^, , (9) 



0( = y and T — g'^^T^j, is the trace of the matter stress-tensor. Note that the parameter 
a is related to inverse of the BD parameter cu. In Einstein frame, the vacuum sector of 
the action consists of a scalar field minimally coupled to Einstein's gravity. The important 
difference between the Einstein frame representation of BD model and minimally coupled 
scalar field models is that in the former the scalar field interacts with matter systems. This 
anomalous gravitational coupling has no counterpart in Einstein's gravity. It implies that the 
stress-tensors of matter and the scalar field are not separately conserved. This can be easily 
checked by applying Bianchi's identities to (7) which leads to 

WT,, = -WT^^, = ^aT V.^ (10) 

The parameter a measures the strength of the interaction. Here there are two important 
points in order : First, the parameter a is positive « > 0. It implies that energy transfer is 
from scalar field ip to matter systems. This feature is consistent with the second law of ther- 
modynamics [16]. Second, the model (4) should be constrained by local gravity experiments 
to avoid violation of weak equivalence principle. It is well-known that these constraints are 
satisfied when a; >> 1 [9]. In the model (4) this translates into a << 1. This means that the 
theory can pass local tests if interaction of the scalar field ip with matter fields is sufficiently 
small. We will return to this issue later. 

We apply the field equations (7) and (8) to a spatially fiat Friedmann-Robertson- Walker space- 
time 

ds'' ^ -de + a\t){dx'' + dy'' + dz'') (11) 
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with a{t) being the scale factor. To do this, we take the matter system to be a pressureless 
perfect fluid (dust) with energy density pm- In this case, the gravitational equations (7) give 

3^ = k{pm + P^) (12) 

2- + — ^-kp^ 13 

where k — SttG, p^ — |<^^ + U{(p) and — |<^^ — U{(p). We may use the first equation to 
rewrite the second one as 

d d^ 1 

a a"^ 2 



+ 2- = -kpm + kU{ip) (14) 



The equation (8) gives 



^d . dU((p) 1 

ip + 3-cp + —j^ ^ --a Pm 15 
a d(fi 2 

On the other hand, the conservation equations (10) become 

d 1 

Pm + 3-pm = if Pm (16) 

p^ + 3-(a;^ + \)p^ = ~a (p pm (17) 

with bj^p = Pip/ Pip being the equation of state parameter of the scalar field (p. The equation 
(16) can be solved which gives the following solution 

Pma^ = Pmoe^"^ (18) 

where pmo is the present matter energy density in the universe. To proceed further, we intro- 
duce an ansatz 

(^=Aina (19) 

in which /3 is a positive constant parameter of order of unity. One of the advantages of this 
ansatz is that it brings the solution (18) into the following form 

Pm = PmOa"'+' (20) 

where e = aP/2y/k — /3[2(2uj + 3)]~^ > 0. This is similar to the rule presented by some 
authors for characterizing decaying law of vacuum energy into dark matter [17]. It states that 
the scalar field (p is constantly decaying into the matter so that the latter will dilute more slowly 
compared to its standard evolution pm oc a~^. Since the observational lower bound imposed 
by solar system experiments on the BD parameter is a; >> 1, we should have £ << 1 which 
means that evolution of matter density has a small deviation with respect to the standard one 
in Einstein's gravity. 

If we put the ansatz (19) into (15), we obtain 

d ^d^ Vk ,dU((p) 1 , 

- + 2^ = -^{^^ + T^aPm) 21 

a p (fi 2 
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Comparing the latter with (14) leads to a consistency relation 

'-^ + VkpUM = -fi\f + .)/>™.e^'-* (22) 

This consistency relation will be satisfied for an appropriate potential U{ip). To find the form 
of the potential, we solve this first order differential equation which gives the following solution 

U{<p) = -^p^oe^^^-'^^ + Ce-^^^^ (23) 



where 7 = ^-ij^^l^ and C is an integration constant. Thus the relation (19) is a solution of the 
field equations for an exponential potential of the form (23). This double exponential potential 
is similar to the potential which is used in some quintessence models [18] [19]. It is shown that 
this kind of potential of the quintessence field can lead to solutions which first enter a period 
of scaling through the radiation and matter domination eras and then smoothly evolve to 
dominate the energy density for a wide range of initial conditions of the field^[18]. Moreover, 
single exponential potentials are popular in modified f{R) gravity models [20]. These models 
are conformally equivalent to BD models with potentials which their forms are closely related 
to the functional form of the f{R) functions [21]. In that context, single exponential potentials 
correspond to power law f{R) gravity models [22]. 

The integration constant C can be determined by noting the fact that when we set = 0o ~ 
in the action (1), then V{(p) characterizes the vacuum energy density corresponding to a 
cosmological constant, namely V{(t)) — A/G. In this case, cp — and then 

U{ip = 0) = V{cl>{^)) = p^o (24) 

with p^Q being the vacuum energy density in the Einstein frame. Applying the latter condition 
to the relation (23) gives C = p^o + ■jpmo- The relation (23) takes then the form 

Ui<p) = Wmoie--^'^^ - ef^^-'^n + p^oe-^^^ (25) 
For this potential function, the Friedmann equation (12) becomes 

3 
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;i - 7)fi„oa-'+^ + {n^o + 7^mo)a-^ ] (26) 



where flmo — Pmo/Pc, ^<fO — Pipo/Pc and pc — SH^/k is the critical density. From the equations 
(13) and (26), it is straightforward to show that the deceleration parameter q — —l — -^ takes 
the following form 

- i{(i + - (3 - - ^_(^^^);,^i)-3..... i-'} 

where we have used a{z) = {z + 1)~^. This relation gives deceleration parameter in terms of 
the redshift and constant free parameters u and /3. We plot evolution of q{z) in fig.l. The 



^'Note that Einstein frame representation of BD models are effectively equivalent to the so-called coupled 
quintessence models in which the quintessence field interacts with matter sector. 
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Figure 1: The plot of deceleration parameter q{z), given by (27), for some values of the 
parameters u and /3. We have set = |. The three lines correspond to e = 0.147, 0.015, 0.002 



from top to bottom. 



figure shows that accelerating expansion is possible for positive small and large values of the 
BD parameter. However, it should be noted that a small u is excluded by two independent 
observations. First, the local gravity tests which require that u > 40000 [9]. Second, from 
the relation (20) we must have e << 1 which corresponds to w >> 1. Otherwise the universe 
will expand accelerated in the matter dominated era, which is against the observation of SNe 
la that our universe expanded decelerated before the redshift z ~ 0.5 [23]. To clarify this 
point, we consider the Friedmann equation (12) or (26) in matter domination regime in which 
Pm >> Pip- In this case, one can simply check that a oc t^/^~^ and H = -^t~^. Only for 
6 << 1, we can expect that the model leads to a decelerating expansion in matter-dominated 
regime whose existence is also fundamental for the structure formation process to take place. 
It is worthwhile to compare our results with recent observations. To do this, we consider a 
parametric approximation of the deceleration parameter along the cosmic evolution, given by 



The analysis is performed using the recent SNe la observational data, the so-called Union2 
sample of 557 events [27]. These observations constrain the parameters go and qi as follows 
[24] [26] : for vanishing spatial curvature, go = —0.66 ± 0.03(lcr) ± 0.07(2cr) and gi = 1.54 ± 
0.19(lcr) ± 0.38(2cr). In fig.2, we show the evolution of the deceleration parameter q{z) for the 
resulting 2cr intervals of go and gi obtained using the data set considered above. As the figure 
indicates, these observations give a bound on the redshift at which the universe switches from 
decelerated to accelerated expansion. We may use these bounds on the transition redshift to 
constrain our model parameters. However, it should be noted that between the two parameters 
u and P appeared in the expression (27) the former has been already constrained by local 
gravity experiments. Thus for a given u, (27) gives the deceleration parameter in terms of 



[24] [25] [26] 
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Figure 2: Evolution of the deceleration parameter q{z) for the resulting 2a intervals of go and 
qi obtained using the data set provided by Union2 sample. 



the redshift and the parameter /3. Exploring the equation (27) reveals that, in this case, for 
a particular value of /3 the transition redshift takes its maximum value. This behavior is 
indicated in fig. 3 for u = 40000. This figure shows that for this value of the BD parameter the 
maximum transition redshift lies within the range given by the Union2 sample. This constrains 
the parameter (3 to be up to 0.7. 

In summary, we have considered the possibility that a self-interacting Brans-Dicke field 
accounts for the accelerated expansion of the Universe in Einstein conformal frame. In our 
analysis, the key point is the ansatz (19) which has two important features. First, it provides 
a class of solutions of the field equations in terms of a potential of the form (23). We have 
shown that these solutions are consistent with late-time accelerating expansion of the universe 
for large values of the BD parameter u. Second, it modifies the evolution of matter density to 
(20) which is the simplest possible way of stating that the matter dilution is attenuated due 
to its interaction with the scalar field ip. This evolution law indicates that the deviation from 
the standard evolution is characterized by a positive constant parameter e which quantifies 
the decay rate. Our analysis also indicates that recent accelerating expansion is possible for 
e « 1 (or a; >> 1) in accord with local gravity tests. 
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Figure 3: The plot of deceleration parameter in terms of /3 and z. 
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